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Transmit Optimization with Improper Gaussian 
Signaling for Interference Channels 

Yong Zeng, Cenk M. Yetis, Erry Gunawan, Yong Liang Guan, and Rui Zhang 

Abstract 

This paper studies the achievable rates of Gaussian interference channels, when the improper or 
circularly asymmetric complex Gaussian signaling is applied. For the Gaussian multiple-input multiple- 
output interference channel (MIMO-IC) with the interference treated as Gaussian noise, we show that the 
user's achievable rate can be expressed as a summation of the rate achievable by the conventional proper 
or circularly symmetric complex Gaussian signaling in terms of the users' transmit covariance matrices, 
and an additional term, which is a function of both the users' covariance and pseudo-covariance matrices. 
The additional degrees of freedom given by the pseudo-covariance matrix, which is conventionally set 
to be zero for the case of proper Gaussian signaling, provide an opportunity to improve the achievable 
rates of Gaussian MIMO-ICs by employing the improper Gaussian signaling. To this end, widely linear 
precoding is proposed, which efficiently maps the proper information-bearing signals to the improper 
transmitted signals for any given pair of transmit covariance and pseudo-covariance matrices. In particular, 
for the case of two-user Gaussian single-input single-output interference channel (SISO-IC), we propose 
a joint covariance and pseudo-covariance optimization algorithm with the improper Gaussian signaling, 
which achieves the Pareto-optimal rates. A suboptimal low-complexity algorithm with separate covariance 
and pseudo-covariance optimizations is also proposed, which guarantees a rate improvement over the 
conventional proper Gaussian signaling. Numerical results are provided to verify the performance of our 
proposed algorithms as compared to existing benchmark schemes. 

Index Terms 

Improper Gaussian signaling, interference channels, pseudo-covariance, widely linear precoding, 
optimization. 

I. Introduction 

The interference channel (IC) models multi-user communication systems where each transmitter is 
intended to send independent information to its corresponding receiver while causing interference to all 
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other receivers. Although information-theoretic study of the IC has a long history [1], characterization of 
its capacity region still remains an open problem in general, except some special cases such as that with 
the presence of "strong" interference |2j3 1. For the single-input single-output IC, termed SISO-IC, the best 
achievable rate region to date is obtained by the celebrated Han-Kobayashi scheme (3j. Recently, it has 
been shown in [4] that a particular form of this scheme achieves within one bit to the capacity region of 
the two-user SISO-IC. Since such capacity-approaching techniques require multi-user encoding/decoding, 
which are difficult to implement in practical systems, a more pragmatic approach is to employ single- 
user encoding and decoding by treating the interference as Gaussian noise at the receivers. In fact, this 
simplified approach has been shown to be sum-capacity optimal when the interference is below a certain 
threshold |5|g| . 

Under the assumption of single-user detection (SUD) with the interference treated as Gaussian noise, the 
transmit optimization problem for Gaussian ICs reduces to resource allocation among the transmitters for 
interference mitigation, which has received significant attention in the last few decades. Early works on re- 
source allocation for ICs mostly focused on power control since single-antenna terminals were considered 
(see e.g. @-||9j and references therein). When transmitters/receivers are equipped with multiple antennas, 
the system performance can be improved via transmit/receive beamforming. One useful technique applied 
for optimizing transmit beamformers is to transform the design problem into an equivalent receiver 



beamformer optimization problem via the so-called uplink-downlink or network duality principle [10|- 



p4| . The transmit beamforming optimization problems for power minimization and signal-to-interference- 
plus-noise ratio (SINR) balancing can also be directly solved by convex optimization techniques, such 
as the second-order cone programming (SOCP) [15] and semidefinite programming (SDP) [16]. 

For Gaussian ICs with the interference treated as Gaussian noise, the sum-rate maximization problem 



is in general difficult to be solved globally optimally due to non-convexity. In |17|, it was shown 
that finding the globally optimal beamformer for weighted sum-rate maximization (WSRMax) in the 
Gaussian multiple-inut single-output IC (MISO-IC) is a NP-hard problem. Algorithms based on the 



principle of interference pricing have been proposed for achieving the local optimum |18|. In |19|, a 
distributed algorithm was proposed for the MISO-IC by using the virtual SINR framework. Gradient 
descent algorithms have been proposed for the Gaussian multiple-input multiple-output IC (MIMO-IC) 
over the transmit co variance matrices p0| or the precoding matrices pT) . More recently, for Gaussian 
SISO-IC, single-input multiple-output IC (SIMO-IC) and MISO-IC, the globally optimal solutions to 
WSRMax problems have been obtained under the monotonic optimization framework p2j-p6*J. However, 
the complexity of such globally optimal algorithms increases exponentially with the number of users, and 
the generalization to the more general MIMO-IC remains unknown. An alternative technique for solving 
the WSRMax problems for MIMO-ICs is via iteratively minimizing the weighted mean-square-error 
(MSE), which utilizes the inherent relationship between the mutual information and MSE | 27|28 1. The 
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study of ICs with game-theoretic models has also been given in 1 29 1 and references therein. Moreover, it 
is worth mentioning that there has been a great deal of interest in the last few years on studying Gaussian 
ICs from the degrees-of-freedom (DoF) perspective [30]. A key technique to achieve higher DoFs than 
previously believed for ICs is known as "interference alignment (IA)" [ 31|32 |. Since DoF only gives 
the capacity approximation at asymptotically high signal-to-noise ratio (SNR), a number of IA-based 
precoding schemes with improved sum-rate performance at practical SNRs have been proposed in (2TJ, 
(33j-(36j. 

Regarding characterization of the achievable rate region for ICs with interference treated as noise, 
various solutions have been obtained for the SISO-IC (37), SIMO-IC [24], and MISO-IC (38)-|j4T). The 
Pareto boundary of the achievable rate region for the IC consists of all the achievable rate-tuples at 
each of which it is impossible to improve one user's rate, without simultaneously decreasing the rate 
of at least one of the other users. Traditional approaches for characterizing Pareto boundaries of ICs 



are via solving WSRMax problems. However, as pointed out in |39|, the WSRMax approach cannot 
guarantee the finding of all Pareto-boundary points due to the non-convexity of the achievable rate set. 
An alternative method based on the concept of rate-profile was thus proposed in [39], which is able to 
characterize the complete Pareto boundary for ICs. Besides, the rate-profile approach generally results in 
optimization problems that are easier to handle than conventional WSRMax problems for ICs [ 24|39 1 . 

It is necessary to point out that in all the aforementioned works on Gaussian ICs, the transmitted signals 
are assumed to be proper or circularly symmetric complex Gaussian (CSCG) distributed. A key property 
of proper Gaussian random vectors (RVs) is that their second-order statistics are completely specified by 
the conventional covariance matrices under the zero-mean assumption. In contrast, for the more general 
improper Gaussian RVs, an extra parameter called pseudo-covariance matrix is required for the complete 
second-order characterization |42|-||44}. Most of existing works on Gaussian ICs have adopted the proper 
Gaussian assumption without any justification. This is probably due to the common practice of modeling 
the additive receiver noise as proper Gaussian, and the well-known maximum-entropy theorem [42], i.e., 
proper Gaussian RVs maximize the differential entropy for any given covariance matrix. As a result, proper 
Gaussian signaling has been shown to be capacity optimal for the Gaussian point-to-point, multiple-access 
(MAC) and broadcast channels (BC). However, for Gaussian ICs with the interference treated as noise, 
the improper Gaussian signaling provides a new opportunity to improve the achievable rates over the 



conventional proper Gaussian signaling [45]. For instance, it was shown in [46 1 that improper Gaussian 
signaling, together with IA, is able to improve the DoF for the three-user SISO-IC with time-invariant 
channel coefficients at the asymptotically high SNR. In | 47|48 1, it was shown that the achievable rate 
region can be enlarged with improper Gaussian signaling even for the two-user SISO-IC at finite SNR. 
This is interesting since it is known that for the two-user SISO-IC, no DoF gain is achievable with IA |45j. 
Specifically, a rank-1 signaling scheme was proposed in (47]], where the transmit covariance matrices for 
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the equivalent real- valued MIMO-ICs are restricted to be rank one. In |48|, the Pareto-optimal transmit 
covariance matrices are obtained by an exhaustive search method. 



The prior works |46|-|48] on the study of improper Gaussian signaling for ICs are based on the 
equivalent double-sized real-valued MIMO-IC matrix by separating the real and imaginary parts of 
each element in the complex-valued MIMO-IC matrix. Although any complex-valued system can be 
transformed into an equivalent real- valued system, as pointed out in [ 43|49 |, much of the elegancy of 
the system description is lost. Therefore, in this paper, we adopt the complex-valued channel model 
for studying the improper Gaussian signaling in Gaussian ICs to gain some new insights. The main 
contributions of this paper are summarized as follows: 

• Based on existing results on improper Gaussian RVs, we derive a new achievable rate expression 
for the general .fT-user MIMO-IC, when the improper Gaussian signaling is applied. Our result 
shows that the user's achievable rate can be expressed as a summation of the rate achievable by the 
conventional proper Gaussian signaling in terms of the users' transmit covariance matrices, and an 
additional term, which is a function of both the users' covariance and pseudo-covariance matrices. 
This new result implies that the use of improper Gaussian signaling for MIMO-ICs with interference 
treated as noise is guaranteed to improve the achievable rate over the conventional proper Gaussian 
signaling with any given set of covariance matrices of transmitted signals, by further optimizing 
their pseudo-covariance matrices. 

• For any given pair of signal covariance and pseudo-covariance matrices at each transmitter, we 
consider the practical problem of generating improper transmitted signals from proper information- 
bearing signals. Based on existing techniques on improper Gaussian signal processing | 50|51 1, we 
propose an efficient method for this implementation, named as widely linear precoding. 

• By adopting the rate-profile method, we formulate the optimization problem for the two-user SISO-IC 
to characterize the Pareto boundary of the achievable rate region with improper Gaussian signal- 
ing. A joint covariance and pseudo-covariance optimization algorithm is proposed, which achieves 
Pareto-optimal rate-pairs. A suboptimal and low-complexity algorithm with separate covariance and 
pseudo-covariance optimizations is also proposed, which guarantees the rate improvement over the 
conventional proper Gaussian signaling with any transmit covariance. 

The rest of this paper is organized as follows. Section |n] studies improper Gaussian signaling for 
the general MIMO-IC, where the achievable rate expression is derived and the widely linear precoding 



is proposed. Section III focuses on the two-user SISO-IC setup, where the problem formulation for 
characterizing the Pareto-boundary of the achievable rate region is given. In Section [iVj joint covariance 
and pseudo-covariance optimization algorithm for the two-user SISO-IC is proposed. In Section [Vj a 
suboptimal but more efficient algorithm with separate covariance and pseudo-covariance optimization is 



presented. Section VI presents numerical results. Finally, we conclude the paper in Section VII 
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Notations: In this paper, scalars are denoted by italic letters. Boldface lower- and upper-case letters 
denote vectors and matrices, respectively. I denotes the identity matrix and denotes an all-zero matrix. 
For a square matrix S, Tr(S), |S|, S _1 denote the trace, determinant and inverse of S, respectively. 
S >z and S >- mean that S is positive semidefinite and positive definite, respectively. C MxN and 
!>Mx7V denote the space of M x N complex and real matrices, respectively. For an arbitrary matrix A, 
A*, A T , A^ and rank{A} represent the complex-conjugate, transpose, conjugate transpose and rank of 
A, respectively, diagjx} represents a diagonal matrix with the elements in the main diagonal given by 
x. blkdiag{ A, B} represents a block diagonal matrix with A and B on the diagonal blocks. ~ means 
"distributed as". CA/"(x, S) represents the CSCG RV with mean x and covariance matrix 5]. The symbol 
i represents the imaginary unit, i.e., i 2 = —1. and represent the real and imaginary parts of 
complex numbers, respectively, /(x; y) represents the mutual information between two RVs, x and y. 

II. Improper Gaussian Signaling for MIMO-IC 

Consider a iT-user MIMO-IC, where each transmitter is intended to send independent information 
to its corresponding receiver, while interfering with the other K — 1 receivers. Denote the number of 
transmitting and receiving antennas for each user by M and N, respectively. Assuming the narrow-band 
transmission, the equivalent baseband signal model at each receiver can be expressed as 

y fc (n) = H fc fcx fc (n) + ^H fcj x 3 (ra) + n fc (n), k = 1, • • • , K, (1) 

where n is the symbol index, £ £_NxM denotes the direct channel matrix from transmitter k 

to receiver k, while Hkj,j / k, denotes the interference channel from transmitter j to receiver k; 
we assume quasi-static fading and thus all channels are constant over re's in ([T]) for the case of our 
interest; rifc(n) represents the independent and identically distributed (i.i.d.) CSCG noise vector with 
nfc(n) ~ CAf(0,a 2 I); and Xfc(ra) £ C M is the transmitted signal vector from transmitter k, which is 
independent of Xj(n) for j / k. In this paper, for the purpose of exposition, we assume that symbol 



extensions over time as in [46] are not used. Hence, Xfc(n) is independent over n. For brevity, n is omitted 
in the rest of this paper. Different from the conventional setup where the proper Gaussian signaling is 
assumed, i.e., x^ ~ CJ\f(0, C Xfc ), in this paper, we consider the more general improper Gaussian signals, 
for which some preliminaries are given next. 

A. Preliminary for Improper Random Vectors 

For a zero-mean RV z £ C n , the covariance matrix C z and pseudo-covariance matrix C z are defined 



as |42 1 



C z = E{zz H ), C z = E(zz T ). (2) 
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By definition, it is easy to verify that the covariance matrix C z is Hermitian and positive semidefinite, 
and the pseudo-covariance matrix C z is symmetric. 



Definition 1. [421: A complex RV z is called proper if its pseudo-covariance matrix C z vanishes to a 



zero matrix; otherwise it is called improper. 

The condition for uncorrelatedness of two complex-valued RVs are specified as follows. 



Lemma 1. [42]: Two zero-mean complex RVs x and z are uncorrelated if and only if C xz = and 
C xz = 0, where C xz = E(xz^) and C xz = E(xz T ). 

A more restrictive definition than properness is known as circularly symmetric, which is defined as 
follows. 



Definition 2. [431: A complex RV z is called circularly symmetric if its distribution is rotationally 



invariant, i.e., if z and z = e ia z have the same distribution for any real value a. 

For a circularly symmetric RV z, we have 

C z = C z = E(zz T ) = e i2a C z , Va, 

which implies C z = 0. Thus, circularity implies properness, but the converse is not true in general. 
However, if z is a zero-mean Gaussian RV, then properness and circularity are equivalent, as given by 
the following lemma. 



Lemma 2. 143]: A complex zero-mean Gaussian RV z is circularly symmetric if and only if it is proper. 



For example, the commonly adopted assumption that the noise vector in ([T} is zero-mean CSCG is 
equivalent to that is a proper Gaussian RV, whose pseudo-covariance matrix satisfies C nfc = 0. For 
an arbitrary complex RV z, define C z as the covariance matrix of the augmented vector [z T (z*) T ] T , 
i.e., 



C z = E 



z 




z 


> 


C z C z 


_z*_ 




z*. 




c* c z . 



(3) 



The augmented covariance matrix C z obviously has some built-in redundancy for the second-order 
characterization of z; however, it is useful as shown in the following two theorems. 

Theorem 1. [43]: C z and C z are a valid pair of covariance and pseudo-covariance matrices, i.e., there 
exists a RV z with covariance and pseudo-covariance matrices given by C z and C z , respectively, if and 
only if the augmented covariance matrix C z is positive semidefinite, i.e., C z ^ 0. 
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Note that the conditions of the covariance matrix C z being Hermitian and positive semidefinite, and 
the pseudo-covariance matrix C z being symmetric are already implied by C z >z 0. Furthermore, for the 
improper complex Gaussian RVs, the differential entropy is in general a function of both the covariance 
and pseudo-covariance matrices, which can be directly expressed in terms of C z as shown by the following 
theorem. 

Theorem 2. [43]: The differential entropy of a complex Gaussian RV z € C ra with augmented covariance 
matrix C z is 

h(z)= 1 -log((7re) 2n \C z \). (4) 
Theorem [2] is a generalization of the entropy formula for proper Gaussian RVs: If C z = 0, (|4]) reduces 



to the well-known expression for proper Gaussian RVs, i.e., h(z) = log((7re) n |C z |) W2f. 



B. Achievable Rate with Improper Gaussian Signaling 

In this subsection, we derive the achievable rate by the improper Gaussian signaling for the K- 
user MIMO-IC defined in ([T]). Denote the covariance and pseudo-covariance matrices of the zero-mean 
Gaussian signal vector Xj by C Xj and C Xj , respectively, i.e., 

C Xj . = E{x jX f }, C Xj = E{x,xJ}. 

Since Xj and xj are independent for I ^ j, then by using Lemma [T] and the fact that independence and 
uncorrelatedness are equivalent for Gaussian RVs, the covariance and pseudo-covariance matrices of the 
received signal vector y&, k = 1 , • • • , K, can be obtained as 

K 

C yk = E{y fc yf } = H fcj -Cx,Hg + a 2 I, (5) 

i=l 

K 

C yk = E{y fc yfl} = H,.,C X (6) 
i=i 

where in Q, we have used the fact that the pseudo-covariance of the CSCG noise vector is a zero 
matrix. It is obvious from ((5]) that C Yk is nonsingular. Then with the augmented covariance matrix C yfc 
defined as in ([3]) and using the Schur complement [52} , we obtain 



I Wfc 1 1 ^y k ~ ^y fc Wi- 1 

| C yJ I 1 ~ Cy^CytCyJCyJ, (7) 
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where for an invertible matrix A, Ar T = (A _1 ) T = (A T ) _1 . In the last equality, we have used the 
identities |A| = |A T |, |AB| = |A||B|, and the facts that C yk is Hermitian and C yfc is symmetric. With 
the input signal Xj being Gaussian, the received signal is also Gaussian. Then based on Theorem [2j 
the differential entropy of y k G is given by 

h(y k ) = log(vre) w |C yfc | + |log|l - C y iC yfc C y Jc£ | . (8) 

Denote as the interference-plus-noise term at receiver k, i.e., 

Sfc = HfcjXj + n k . 

Then the covariance and pseudo-covariance matrices of are given by 

C Sfc =^H fcj C Xj Hg +C r 2 I, (9) 

C Sfc =]TH fcj C Xj H^, (10) 

Similarly as for y k , the differential entropy of s k can be obtained as 

h(s k ) = log(7re)"|C. fc | + ilog|l - C- 1 C Bfc C- r C^|. (11) 

Under the assumption that interference is treated as Gaussian noise, the achievable rate at receiver k with 
improper Gaussian signaling can be obtained as 



R k =I(x fc ;y fc ) = h(y k ) - h(y k \x k ) = h(y k ) - h(s k ) 
\<r 2 l + J2f=i H fciC Xj Hg| 



log 



— Rk, proper ({Cxj }) 

+ -wJ ^ yfc y* y*i (12) 

The above equation shows that with improper Gaussian signaling, the achievable rate can be expressed 
as a summation of two terms. The first term, denoted by i?fc, pr oper({C x .}), is the rate achievable by 
the conventional proper Gaussian signaling, which is a function of the transmit covariance matrices 
only. The second term, is a function of both the transmit covariance and pseudo-covariance matrices. By 
setting C Xj = 0, Vj, the second term vanishes and ( fT2] l reduces to the rate expression for the conventional 
case of proper Gaussian signaling. The separability of the achievable rate by improper Gaussian signaling 
provides a general method to improve the achievable rate over the conventional proper Gaussian signaling, 
i.e., for any given covariance matrices obtained by existing proper Gaussian signaling schemes, the rate 
can be improved with improper Gaussian signaling by choosing the pseudo-covariance matrices that 
make the second term in ( fT2] ) strictly positive. It is worth noting that this property does not exist if we 
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convert the complex-valued system in ([T} to an equivalent real-valued system by doubling the input/output 
dimensions. 

In this paper, we are interested in characterizing the achievable rate region with improper Gaussian 
signaling. The achievable rate region for the K-user MIMO-IC consists of all rate-tuples for all users 
that can be simultaneously achieved under a given set of transmit power constraints for each transmitter, 
denoted by Pj,j = 1, ...K, i.e.,: 



R= |J |(ri,-" ,r K ) :0<r k < R k ,k = 1, ...k\, 



(13) 



Tr{C Xj }<P 3 



where R k is given in ( [12] ); C x is the augmented covariance matrix defined in ([3]). The constraint C x >z 
follows from Theorem [TJ 

Definition 3. ft38y: A rate-tuple (r%, ■ ■ ■ , tr) is Pare to optimal if there is no other rate-tuple (r^, • • • , r' K ) 
with (r[, ■ ■ ■ , r' K ) > (n, • • • , Tk ) (r^, • • • , r^) 7^ (n, • • ■ , r^-), where the inequality is component- 
wise. 



In Section III IV] we will consider the transmit covariance and pseudo-covariance optimization for 



achieving the Pareto rate boundary for the special case of two-user SISO-IC. 

C. Widely Linear Precoding 

In this subsection, we consider the practical problem of how to efficiently generate the transmitted 
signal at each transmitter given any pair of covariance matrix C Xj and pseudo-covariance matrix C X] , 
from an information-bearing signal dj that is selected from CSCG (proper Gaussian) codebooks. Without 
loss of generality, assume that d, ~ CAf (0,1), i.e., the covariance and pseudo-covariance matrices of dj 
are given by 

C d] =1, C dj =0, j = !,■■■ ,K. (14) 

First, consider the conventional linear precoding for a multi-antenna transmitter with proper Gaussian 
signaling, 

Xj = Vjdj, (15) 

where Uj is the precoding matrix. Then the pseudo-covariance matrix of Xj is given by C Xj = Uj U J = 
0. This implies that the conventional linear precoding is not able to map the proper Gaussian signal dj 
to the improper Gaussian signal Xj. 

Since the augmented covariance matrix defined in ([3]) contains both the covariance and pseudo- 
covariance matrices, a sufficient condition for a RV Zj to have covariance matrix C Xj and pseudo- 
covariance matrices C Xj is that its augmented covariance matrix satisfies C Zj = C Xj . This is ensured by 
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the following transformation 



(16) 



where C£ 3 denotes the generalized Cholesky factor of the positive semidefinite matrix C x ., which is 
defined by C x = C| (Cl,)^ [441. Since C d = I due to ([14]), it is easy to verify that Zj in ([16]) satisfies 

3 3^ 

C z = C x .. A common method for finding Cj is via eigenvalue decomposition (EVD). Specifically, let 

3 3 J ^ ^ 

the EVD of C x be expressed as C x = UDU H ; then C Yj = UDs is obtained. However, it is worth 

3 3 ^ 

pointing out that the above obtained C^L cannot satisfy ( [To] ) in general. This is because the two vectors 
zj and z* in ( fTB"] ) are complex conjugate of each other; therefore, the transformation matrix C x should 
be designed with more care as compared to the conventional EVD. Furthermore, if we can find an C Xj 
such that it has the following structure: 



Bi 



B 2 



then it follows that (16) is equivalent to the following two sets of equations 



Bjd, +B 2 d*, 



z* = B* 2 dj +Bld*. 



(17) 

(18) 
(19) 



It is easy to verify that the two equations given in ( [T8] ) and ( [T9| ) are consistent, i.e., ( fT9] ) is simply obtained 



by taking the complex conjugate on both sides of ( [18] ) and vice versa. Therefore, the remaining task is 
to find one Cj^ with the structure given by (T7\ . To achieve this end, we define the following 2M x 2M 
unitary matrix ||49j: 

_ A 1 



V2 



a ' 

-il 



, TT 



H 



I. 



For any real-valued matrix A G 



p2Mx2M 



, it can be easily verified that the matrix TAT^ has the specific 



structure given in (17), i.e., 



TAT 



H 



A* A* 



Theorem 3. [43 491 The augmented EVD of the augmented covariance matrix C x . £ C 
in ([3]) is given by 



2Mx2M 



V A V 



H 



(20) 



defined 



(21) 



with 



V = TVT ff , A = TAT H , 



July 24, 2012 



DRAFT 



11 



where V S l^** 2 ™ is a real-valued orthogonal matrix, A = blkdiag{Ai, A2}, Ai = diag{Ai, A3, • • • , A2M— 1} 
arcc? A2 = diag{A2, A4, • • • ,A2m} consisting of the eigenvalues of C Xj w/f/i f/ze ordering X\ > A2 > 

• • • > A2M- 



For any given C Xj ^ 0, all of its eigenvalues are nonnegative, i.e., Aj > 0, 1 = 1, 
can be written as 



,2M. Then (21 



C x = TVAV a T H = (TVA 1/2 T a )(TVA 1/2 T H ) ff . 



Thus we have 



Cl = T(VA 1/2 )T fl 



B 2 
B^ 



(22) 



where the last equality follows from (20) and the fact that VA 1 / 2 is a real- valued matrix. 

From ( [16] ) and ( [22] ), it follows that to obtain the transmitted signal vector Xj, which is generally 
improper with the covariance and pseudo-covariance matrices specified by C Xj , the following precoding 
needs to be applied to the proper information-bearing signal dj\ 

Xj - = B 1 d i + B 2 d*, (23) 

where Bi and B2 are the corresponding blocks of TVA 1 / 2 T^ as shown in (|22]). Following similar 



terminologies used in existing literatures on improper signal processing such as [50], [51 1, we refer to 



the precoding given in ( [23] ) as widely linear precoding. Note that if B2 = 0, which is the case when C x 
is block-diagonal (i.e., C Xj = 0), ( |23[ ) reduces to the conventional linear precoding for proper Gaussian 
signaling given by ( [T5] ). Lastly, in terms of the real-valued representation, ( |23] ) can be re-expressed as 







-9{xi}- 





5t{Bi+B 2 } 9{B 2 -Bi} 
9{B 2 + Bi} K{Bi-B 2 } 



R{d,-} 



III. Pareto Boundary Characterization for the Two-user SISO-IC 

The Pareto boundary of the achievable rate region consists of all the Pareto-optimal rate-tuples as 
shown in Definition [3] For the conventional setup where all the transmitted signals are restricted to be 
proper Gaussian, i.e., C Xj = 0, Vj, the Pareto boundary characterization has been obtained for the SISO- 
IC (37|, SIMO-IC (24), and MISO-IC (38j-||4lJ. However, for the MIMO-IC case, the proper Gaussian 
signaling scheme that achieves all Pareto optimal rate-tuples remains unknown. In the remaining part of 
this paper, we focus on the particular two-user SISO-IC, with the objective of characterizing its Pareto 
boundary with improper Gaussian signaling by optimizing both the covariances and pseudo-covariances 
of transmitted signals. The input-output relationship for the two-user SISO-IC can be simplified from ([JJ) 
as 

yi = huxi + h 12 x 2 + n 1 , 



y 2 = h 2 ixi + h 22 x 2 + n 2 , 



(24) 
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where h k j = \h k j\e l ^ kl ,k, j E {1,2}, is the complex scalar channel from transmitter j to receiver /cQ 
Denote the covariance and pseudo-covariance of the transmitted signals 

C Xl = E( Xl xl), C Xl = E(sixi), (25) 
C X2 = E(x 2 x* 2 ), C X2 = E(x 2 x 2 ). (26) 

Note that in the above, C Xl and C X2 are nonnegative real numbers equal to the transmit power values 
of the corresponding users, while C Xl and C X2 are complex numbers in general. From Theorem [TJ it is 
easy to verify that the following conditions are both necessary and sufficient for C Xl and C Xl ( C X2 and 
C X2 ) to be a valid pair of covariance and pseudo-covariance for x\ (x 2 ): 

\C Xl \ < C Xl , \C X2 \ < C X2 . (27) 

The covariance and pseudo-covariance of the received signal y k can be computed as 

C yk = E(y k y* k ) = \h kl \ 2 C Xl + \h k2 \ 2 C X2 + a 2 , (28) 

C yk = E(y k y k ) = h 2 kl C Xl + h 2 k2 C X2 , (29) 
For the interference-plus-noise term s k = h^xj. + n k , where k = {1, 2}\k, we have 

Cs k = \h k - k \ 2 C x ^ + a 2 , C Sk =h 2 k - k C x ^. (30) 



Following similar derivations as Section II-B the achievable rates for the two-user SISO-IC can be 
obtained as 



Rf so = I(x k ;y k ) = h(y k ) - h(s k ) 

i r 2 - \r i 2 

= hog y * l J kl (31) 

= + ) (32) 

V v ' 

A oSlSO \ 
-""fe, proper \° ^1 > 0a: 2 ) 

To characterize the Pareto boundary of the achievable rate region, we adopt the rate -profile technique 



proposed in |39|. Specifically, any Pareto-optimal rate-pair can be obtained by solving the following 



'Since a phase rotation can be applied at each of the receivers with coherent demodulation, without loss of generality, the 
direct channel gains h\\ and /122 can be assumed to be real values. However, this assumption will not change the results of this 
paper. 

2 Note that we still use the terms covariance and pseudo-covariance for scalar random variables, although the more commonly 
used terms are variance and pseudo-variance, respectively. 
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optimization problem with a given rate-profile specified by a= (0:1,0:2): 

(PI): max. _ R 

Cxi )@Xx 

1 c 2 —\C I 2 

s.t. -log-^ -J^- > a k R, k = 1,2, (33) 

2 ^C 2 -\C Sk \ 2 " 

0< C Xj <Pj, i = 1,2, (34) 
\C Xj \<C X] , j = l,2, (35) 

where o^ denotes the target ratio between user k's achievable rate and the users' sum-rate, R. Without 
loss of generality, we assume that > 0, k = 1,2 and a\ + «2 = 1- Denote the optimal solution to (PI) 
as R*, then the rate-pair (a\R*, 02-R*) must be on the Pareto boundary corresponding to the rate-profile 
given by (0.1,0.2). Thereby, by solving (PI) with different rate -profile parameters (01,02), the complete 
Pareto boundary for the achievable rate region can be found. 

IV. Joint Covariance and Pseudo-Covariance Optimization 

In this section, (PI) is solved by jointly optimizing the covariance and pseudo-co variance. First, note 
that if a given target sum-rate R corresponding to the rate-profile (cui, 02) is achievable, then the optimal 
solution to (PI) satisfies R* > R; otherwise, R* < R. Therefore, (PI) can be solved via solving a 
sequence of feasibility problems, each for a given target sum-rate R. Furthermore, R* can be found 



with the simple bisection search method |52|. The remaining task is thus to solve the feasibility problem 



resulting from (PI) for a given R, which can be formulated as 

(Pl.l): Find C Xl ,C X2 ,C Xl ,C X2 

s.t. \C y f + e 2a " R Cl < C 2 k + e 2a * R \C Sk \ 2 yk (36) 

< C Xj < Pj, Vj, (37) 

\C Xj \<C x „Vj. (38) 

Lemma 3. Without loss of generality, we can assume that the pseudo-covariance C Xl is a nonnegative 
real number. 



Proof: By substituting the expressions of C Vk and C Sk into ( |36| ) with p9| ) and ( |30| ), respectively, it is 
easy to see that a common phase rotation to both C Xl and C X2 does not change the feasibility of (Pl.l). 
In other words, if {C Xl ,C X2 ,C Xl ,C X2 } is a solution to (Pl.l), then so is {C Xl , C X2 , g iuj C Xi , c tlJJ C X2 \ for 
any ui. Therefore, without loss of generality, we may choose uj so that C Xl is real and nonnegative. ■ 
Next, we show that (Pl.l) can be transformed into a sequence of SOCP problems by approximating 



the right hand side (RHS) of the constraint in (36). First, a complex variable z k is introduced, with C, 
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and e akR \C Sk \ as the real and imaginary parts, respectively, i.e., 

z k ±C yk +ie a » R \C Sh \, (39) 
Then the RHS of d36b equals to \z k \ 2 . The magnitude of a complex number can be represented as follows: 



Lemma 4. [53] For any complex number z = u + iv, we have 



\z\ = yV + v 2 = max mze'^}, (40) 
/?e[o,i] 

where the maximum is attained when 2tt/3 = arg z, where arg z denotes the phase of z. 

The above result can be easily verified by writing the complex number z in polar form, i.e, z = 
\z\e lSLVgz . An approximation to ( |40] > can be made by evaluating the maximum over a discrete set of 
values for f3, e.g. {3\ = 1/2L, where I € {0, 1, • • • 2L — 1} and L is a positive integer that determines the 
accuracy of the approximation. Specifically, let 

p L (z) ± max mze~ i2 ^ l ' 2L } 
ie{0,-,2L-l} 

= max \Mze- i7Tl / L }\ 
ie{o,- ,L-iy 1 

= max \mz}cos(ln/L) + %{z}sm(lTr/L)\, (41) 

2e{0,-,L-l} 1 

then \z\ is bounded by pl(z) < \z\ < pl{z) sec(ir/2L) |53], where sec(-) denotes the secant function. 
Therefore, by choosing L to be large enough, \z\ can be approximated by pl(z) with a good accuracy. 
In this paper, we choose L = 8, which corresponds an error bound sec(7r/2L) = 1.0196, i.e., the 
approximation error in the worst case is less than two percent. Applying this result to (Pl.l), we can 
approximate the RHS of ([36]) by 

\zk\ ~ max C v , cosU-k/L) + e akR \C s , I sin(Z7r/L) . 

CKKL-l y 11 

The intervals for / can be further narrowed by utilizing the prior knowledge of argzfc. Specifically, it can 
be verified that \C Sk \ < C Sk < C Vk ; thus with z^ defined in (39l, we have argz^ < tan -1 \e akR ) = 



Therefore, / can be restricted to the values such that lir / L < r]k is satisfied for user k, i.e., Ik < rjkL/ir. 
Let 

L k ^\ Vk L/Tr],k = l,2. (42) 

We then have 

~ max Cy k cos(ifc7r/L) + e akR \C Sk | sm(l k Tr/L) . 

0<lk<L k 
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Then for any pair of l\ £ {0, 1, • • • , L{\ and l 2 € {0, 1, • • • , L 2 }, (Pl-1) can be rewritten as 

(P1.2): Find C Xl ,C X2 ,C Xl ,C X2 
s.t. \C y f + e 2a * R C 2 k < 

(C Vk cos(/ fe 7r/L) + e akR \C Sk I sm(l k TT/L)) 2 , Vk 
0<C Xj <Pj, j = l,2, 

\C Xj I < Cj;. , j = 1,2. 



(43) 
(44) 
(45) 



Since the RHS of the constraint in ( |43| ) is a lower bound on that in p6| ), the feasibility of (Pl.l) can be 
obtained as follows: (Pl.l) is feasible if for at least one pair of l 2 }, (PI. 2) is feasible. This statement 
is exact regardless of the value of L. On the other hand, if (PI. 2) is infeasible for all pairs of {li,l 2 }, 
then we claim that (Pl.l) is infeasible. This statement becomes exact when L — > 00. 

In order to solve (PI. 2), denote the magnitude and phase of C X2 by t and 9, respectively, i.e., C X2 = te t6 '. 
Then with Lemma [3] and by substituting C Vk and C Sk with ( |29] > and ( [30] ), respectively, (PI. 2) can be 
recast as 

(PI. 3): Find C Xl ,C x „C Xl ,t,0 



s.t. 



e a * R C S2 



< C Vl cos(lm/L) 

+ e aiR \h 12 \ 2 tsm{hir/L), 

< C y2 cos(/ 2 vr/L) 

+ e a2R \h 2 i\ 2 C Xl sm(l 2 Tr/L), 



(46) 



(47) 



< C Xl < Py, < C X2 < P 2 , 



< C Xl < C Xl 



0<t<C r 



where C Vk , C Sk are affine functions of C Xl and C X2 given in ( |28] ) and ( [30] ), respectively. For any fixed 
value of 6, (PI. 3) is a SOCP problem, which can be efficiently solved by existing software [54 1. Therefore, 
(PI. 3) can be solved by searching over all possible values of 9 until a feasible solution is obtained (or 



not obtained at all). Intuitively, the optimal 9 should minimize the LHS of (J46]» or ( |47] ) to make (PI. 3) 
feasible with the largest margin. Therefore, two good initial values for 9 to test the feasibility of (PI. 3) 
are given by 



9 £ {vr + 2(0ii - 012), vr + 2(<fo 1 - <j) 22 )}, 



(48) 
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where the first value corresponds to the two terms h\ l C Xl and h^te" 10 in (46i being added destructively 



(recall that C Xl is real and cf)kj is the phase of hkj), and the second value makes the two corresponding 
terms in ( |47| ) being added destructively. 

The proposed joint covariance and pseudo-covariance optimization algorithm for the two-user SISO-IC 
is summarized as follows: 



Algorithm 1 Joint Covariance and Pseudo-covariance Optimization for (PI) 
1: Initialize L. Let 0o = {^jplm = 0, • • • D — 1}, where D is the number of evenly spaced 9 values 

to be searched in the interval [0, 2n). Let = {it + 2(<pn — (pw), vr + 2(4>2i — $22)} U 0o- 
2: Let Rl = 0, Ru be a value large enough so that the optimal solution to (PI) satisfies R* < Rjj. 
3: repeat 

4: Set R^ (R L + Ru)/2. Let L k = \L tan~ 1 (e° lkR ) /tt] , k = 1,2. 

5: For 9 € G, l x e {0, • •■ l 2 G {0, • • • ,L 2 }, solve the SOCP feasibility problem (P1.3). If 

feasible, set Rl R and go to step 4; otherwise, repeat step 5. 
6: If (PL 3) is infeasible for all sets of {9, l±, I2}, set Ru <— R. 
7: until Ru — Rl < e., where e is the error tolerance. 



Remark 1. Algorithm^gives the globally optimal solution to (PI) with L — > 00 and D — > 00. However, 
for even moderate values such as L = 8 and D = 10, it provides the solution very close to the optimal, 



as verified by numerical results in Section VI 



V. Separate Covariance and Pseudo-Covariance Optimization 

The algorithm proposed in the previous section, provides the optimal solution to (PI) asymptotically, 
but incurs a very large complexity if high accuracy is required. In this section, by exploiting the fact that 
the user's achievable rate can be separated into two terms in ([32]), we propose a suboptimal algorithm 
with lower complexity for solving (PI) with separate covariance and pseudo-covariance optimizations. 
Specifically, the covariance of the transmitted signals is first optimized by setting the pseudo-covariance 
to be zero, i.e., the proper Gaussian signaling is applied. Then, the pseudo-covariance is optimized with 
the covariance to be fixed as the previously obtained value. With such a separation approach, the obtained 
improper signaling scheme is guaranteed to improve the rate over the optimal proper signaling scheme. 
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A. Covariance Optimization 



When restricted to proper Gaussian signaling with C Xl = and C X2 = 0, by substituting ( |32| ) into 
([33), (PI) reduces to 

(PI. 4): max. r 

s.t. log(l + -^^—) > air, (49) 

0<C Xl < P u < C*, 2 < F 2 . 
For any fixed value r, (PI. 4) can be transformed to the following feasibility problem: 

(PI. 5): Find C Xl ,C X2 
s.t \hn\ 2 C Xl > (a 2 + \h l2 \ 2 C X2 ){e^ r - 1), 
\h 22 \ 2 C X2 > (a 2 + \h 21 \ 2 C Xl )(e a * r - 1), 
< C Xl < Pi, 0<C X2 < P 2 . 

(PI. 5) is a linear programming (LP) problem, which can be efficiently solved by existing techniques such 
as the simplex method or the iterative algorithm in p4| . If r is feasible to (PI. 5), then it follows that 
the optimal solution to (PI. 4) satisfies r* > r; otherwise, r* < r. Thus, (PI. 4) can be efficiently solved 
by solving (PI. 5) together with the bisection method for updating r. 

B. Pseudo-Covariance Optimization 

Denote the optimal solution to the covariance optimization problem (PI. 4) as {r*, C* , C* 2 }. By fixing 
the covariances to C* and C* 2 , (PI) is then optimized over the pseudo-covariances C Xl and C X2 . By 
substituting the first term in the rate expression ([32]) with a^r*, the problem for pseudo-covariance 
optimization is formulated as 

(PI. 6): max. R 

Cxi iC X 2 

s.t. a x r* + -log y \ J 1 >a x R, 

2 l-C s - 2 |C Sl | 2 

a 2 r* + hog — > a 2 R, 
2 1-C7 2 \C S2 \ 2 

Ci 2 /~i*2 \f~i 1 2 /~i*2 
, III i u xi> \^x 2 \ ^ ^ X2 , 
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where C yi , C Sl , C y2 and C S2 are the corresponding covariance terms with the input covariances C* 1 
and C* 2 . Again, if a given R is achievable for certain C Xl ,C X2 , then the optimal solution to (PI. 6) 
satisfies R* > R; otherwise, R* < R. Therefore, (PI. 6) can be solved via solving a set of feasibility 
problems together with the bisection method. It can be easily observed that {C Xl = 0, C X2 = 0, R = r*} 
is feasible to (PI. 6). Therefore, R* > r* is satisfied, i.e., with our proposed separate covariance and 
pseudo-covariance optimizations, the users' sum-rate corresponding to the rate-profile given by {0.1,0.2) 
with improper Gaussian signaling is guaranteed to be no smaller than that obtained with the optimal proper 
Gaussian signaling obtained by solving (PI. 4). Next, we present the algorithm for solving the feasibility 



problem in (PI. 6) for a given R. Substituting ( |29| ) and p0[ ) into (PI. 6) and after some manipulations, 
the feasibility problem for a given R can be formulated as 



(P1.7) min. 

c xi ,C X2 

s.t. ax\h 2 n C Xl + h 2 2 C X2 \ 2 + b 1 < \C X2 \ 2 , 
ci2\h 2 > 1 C Xl + h 22 C X2 \ 2 + 62 < |Cxi| 2 ) 
\C Xl \ 2 <C* 2 



where 



1^X2 I 



(i-i/A)C a 2 



< a 



*2 

x 2 ' 



(51) 
(52) 
(53) 
(54) 



K- k \ 4 



,p k = (&*h(R-r*) t k = l,2,k = {l,2}\k. Since the optimal 



solution to (PI. 6) satisfies R* > r* , we can assume that R > r* in (PI. 7) without loss of optimality. 
Then it follows that j3k > l,&fe > 0,V/c. Next, we show that (PI. 7) can be efficiently solved via solving 
a finite number of SOCP problems. With Lemma [3] we can assume that C Xl is real and nonnegative 
without loss of generality. Again, denote the magnitude and phase of C X2 by t and 6, respectively, i.e., 
C X2 = te ie . Then for any fixed value of 9, (PI. 7) can be transformed into a SOCP feasibility problem 
given by 



(PI. 8): min. 



s.t. 







m.{h 2 n C Xl + h{ 2 te lU ) 
->{h 2 21 C Xl + h 2 22 te id ) 



< t. 



<C X 



c Xl <c* Xi , o<t<a 



x 2 ■ 



Theorem 4. The feasibility problem (PI. 7) can be optimally solved by solving a finite number of SOCP 
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problems (PI. 8), each for a fixed value 9, where 9 can be restricted to the following discrete set: 

9 6 {tt + 2((/> n - <P 12 ),TT + 2(021 - fe)} U 6^ U @ B , 

where 0^4 and are the solution sets for 9 to the following equations: 

[aAhl.Ct +h\ 2 te ie \ 2 + h = t 2 
Q A : { 1 (55) 

(a 2 \h 2 21 C* i +h 2 22 te w \ 2 + b 2 = Cf i 
axl^Xi + h 2 12 C* 2 e i0 \ 2 + h = C* 2 



0g : { " * a ' ^ (56) 

l 22Cx 2 ( 



021^X1 + /ii 2 C*e ifl | 2 + 62 = X? 



Proof: Please refer to Appendix [A] ■ 
Theorem |4]can be intuitively interpreted as follows. For the feasibility problem (P1.7), if the constraint 
( [51) is more "restrictive" than ([52]), then 9 should have a value such that the LHS of ( [51) is minimized. 
This corresponds to 9 = 7r + 2(0n — 0i2) so that h\ l C Xl and h\ 2 C X2 are antiphase. Similar interpretation 
for 9 = it + 2(02i — 022) can be made. On the other hand, if both ( |5"T| ) and ([52]) are equally "restrictive", 



a feasible solution tends to make both constraints satisfied with equality, as given by ( [551 ) an d ( p6| ). 
0_4 and 0g correspond to the cases where either the constraint ( |53| ) or ( |54[ ) is active, which can be 
assumed without loss of generality as shown by Proposition [T] in Appendix [A] The elements in 6.4 and 
Gg can be readily obtained by following the steps in Appendix [B] The proposed separate covariance and 
pseudo-covariance optimization algorithm for (PI) is summarized as follows: 



Algorithm 2 Separate Covariance and Pseudo-covariance Optimization for (PI) 
1: Solve (PI. 4), and denote the optimal solution as {r* ,C* ,C* }. 
2: Initialize Rl = r*, Rjj to the value given by ( [69] ) in Appendix [C] 
3: repeat 

4: Set R <— (Rl+Rjj) /2 and solve the feasibility problem (PI. 7). If feasible, set Rl <— R; otherwise, 
Ru <- R. 

5: until Ru — Rl < e, where e is the error tolerance. 



VI. Numerical Results 

In this section, we evaluate the performance of the proposed algorithms for the two-user SISO-IC 
with numerical examples. Both transmitters are assumed to have the same power constraint P, i.e., 
P 1 = P 2 = p. SNR is defined as P/a 2 . For Algorithm [j] L = 8 and D = 10 are used. 
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Fig. 1: Achievable rate region for the two-user SISO-IC with channel realization and SNR = dB. 
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Fig. 2: Achievable rate region for the two-user SISO-IC with channel realization H^, and SNR = 10 dB. 



A. Rate Region Comparison 

In Fig. [T] and Fig. [2j the achievable rate regions for an example two-user SISO-IC are plotted for 
SNR=0 dB and 10 dB, respectively. The channel matrix for both plots is given by ~HS ' = = 

h,2l h,22 

The proposed improper Gaussian signaling schemes with 



1.5718 - 1.2863* 



■1.2984 + 0.7032* 

L -0.2847 + 0.6700i 0.7802 - 0.6151s 
joint (Algorithm [TJ and separate (Algorithm [2]) covariance and pseudo-covariance optimizations are 

compared with other existing schemes, including the optimal proper Gaussian signaling scheme, the 

optimal improper Gaussian signaling obtained by the exhaustive search method |48| , and the rank- 



July 24, 2012 



DRAFT 



21 



1 improper Gaussian signaling scheme |47|. Both figures reveal that for the given channel H^, the 
achievable rate regions have been significantly enlarged with improper Gaussian signaling over the 
conventional proper Gaussian signaling. The plots also demonstrate that the joint covariance and pseudo- 
covariance optimization algorithm gives almost the optimal solution with only moderate values for L and 
D. Moreover, it is observed that the separate covariance and pseudo-covariance optimization algorithm 
performs close to the optimal as well. It is worth to mention that, even with time-sharing (TS)0 the 
improper Gaussian signaling still outperforms the proper Gaussian signaling, as shown by the dashed 
lines in the two figures. For this particular channel realization, the Pareto boundary points of the achievable 
rate region with TS using improper Gaussian signaling can be obtained by the TS between each of the 
two single-user maximum rate points, and the largest rate corner point by the existing rank-1 scheme 



. The achiev- 



1 47 1 . However, this is not always the case, as illustrated by the next example. 

r-0.8034 + 3.9185i -0.0931 + 0.8952i 

Consider that the channel is given by ~HS > = 

L 0.6488 + 0.4680? -0.7646 + 1.2905i J 
able rate regions for SNR=0 dB and 10 dB are plotted in Fig. [3] and Fig. |4j respectively. It is observed 

from Fig. [3] that for this particular channel realization at SNR=0 dB, there is no performance gain by using 

improper Gaussian signaling over proper signaling. It is also observed that the rank- 1 improper signaling 



scheme |47|, which is based on the equivalent real-valued 2x2 MIMO-IC, gives strictly smaller rate 
region than that by the proper Gaussian signaling. In contrast, our proposed improper signaling schemes 
with either joint or separate covariance and pseudo-covariance optimizations, are observed to perform no 
worse than the optimal proper Gaussian signaling scheme, in accordance with our previous discussions. 
Similar observations can be made in Fig. [4j 

B. Max-Min Rate Comparison 

The rate -profile technique used in characterizing the Pareto boundary of the achievable rate region can 
be directly applied for maximizing the minimum (max-min) rate of the two users without TS. Specifically, 
the max-min problem for the two-user SISO-IC is equivalent to solving (PI) by using the rate-profile 
ol = (1/2, 1/2). An alternative rank-1 max-min solution has also been proposed in [48 1. In Fig. [5j the 



average max-min rates without TS over 100 random channel realizations are plotted for the two-user 
SISO-IC. All the channel coefficients are drawn from the CSCG random variables with zero-mean and 
unit-variance. It is observed that the average max-min rates given by the improper Gaussian signaling 
scheme are significantly higher than that with the proper Gaussian signaling, and the performance gains 
are more pronounced in the moderate to high SNR regime. 

3 The achievable rate region with TS is obtained by taking the convex-hull operation over all the achievable rate-pairs given 
in [TJl. 
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Fig. 3: Achievable rate region for the two-user SISO-IC with channel realization H^ 2 \ and SNR = dB. 
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Fig. 4: Achievable rate region for the two-user SISO-IC with channel realization H( 2 ), and SNR = 10 dB. 



C. Sum-Rate Comparison 

In this subsection, the sum-rate maximization with improper Gaussian signaling is considered. By using 
the equivalent real-valued MIMO-IC of the complex-valued SISO-IC, existing sum-rate maximization 
algorithms in the literature, such as the one via the iterative weighted MSE minimization (WMMSE) 



1 28 1, can be applied directly for maximizing the sum-rate of the two-user SISO-IC when the improper 
Gaussian signaling is employed. However, although the WMMSE algorithm is guaranteed to converge 
to a local maximum of the sum-rate, it is not guaranteed to achieve the global sum-rate maximum. With 
the algorithms proposed in this paper via covariance and pseudo-covariance optimization, we illustrate 
with the following example that our proposed algorithms strictly improve the achievable sum-rate over 
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Fig. 5: Average max-min rate for the two-user SISO-IC. 



that by the WMMSE algorithm. 

In order to apply the WMMSE algorithm to the sum-rate maximization problem when the improper 
Gaussian signaling is applied, we transform the complex-valued channel to the equivalent real-valued 



MIMO channel, similarly as in |46j-||48j. As an example, the iterative WMMSE algorithm in (28J is 
applied to the two-user SISO-IC in (24 1. Denote Qk as the transmit covariance matrix of user k in the 
equivalent real-valued 2x2 MIMO-IC. Without loss of generality, denote the rate-pair obtained by the 
iterative WMMSE scheme by r ■ (61,82), where 81 > 0,8 2 > 0, 81 + 82 = 1. With the rate -profile cx= 
(81, 82), (PI) is solved to obtain a new sum-rate R. If R > r, then the rate -pair obtained by the WMMSE 
algorithm cannot be sum-rate optimal, and a strictly improved sum-rate has been obtained. Table [II] shows 
this result. The equivalence of the channel matrices and the converged transmitter parameters between 
the original complex- valued SISO-IC and the equivalent real- valued MIMO-IC for this example is shown 
in Table ^ with the SNR is as 10 dB. 

VII. Conclusion 

This paper studied the transmitter linear precoding for Gaussian ICs when the improper or circularly 
asymmetric complex Gaussian signaling is applied. Under the assumption that the interference is treated as 
additive Gaussian noise, it was shown that the conventional use of proper or circularly symmetric complex 
Gaussian signaling may cause unnecessary rate losses. A new achievable rate expression for the general 
MIMO-IC was derived, which is expressed as a summation of the rate achievable with the conventional 
proper Gaussian signaling, and an additional term due to the use of improper Gaussian signaling. This 
result provides a useful method to improve the rate over the conventional proper Gaussian signaling by 
separately optimizing the covariance and pseudo-covariance matrices. We also proposed the method of 
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TABLE I: Equivalence between complex-valued SISO-IC and real-valued MIMO-IC 





Complex-valued SISO-IC 


Real- valued MIMO-IC 


Channel 


/in = 2.7388 - 0.2498i, h 12 = 0.9956 + 1.8047i 
h 21 = 0.6680 - 1.6470i, h 22 = 0.4760 + 1.2 706i 


Hn = 
H21 = 


" 2.7388 0.2498" 
-0.2498 2.7388. 
" 0.6680 1.6470" 
-1.6470 0.6680. 


, H12 
■ H22 = 


'0.9956 -1.8047" 
.1.8047 0.9956 _ 
'0.4760 -1.2706" 
.1.2706 0.4760 _ 




WMMSE Init. 


C X1 = 10.000, C X1 = 9.546e i0 5512 
C X2 = 10.000, C X2 = 7.6118e l2 - 8995 


Qi = 


9.0660 2.4998" 
2.4998 0.9340. 


,Q 2 = 


'1.3051 0.91251 
.0.9125 8.6949J 


WMMSE Sol. 


C X1 = 10.000, C X1 = I0e l0 - 6028 
C X2 = 10.000, C X2 = lOe"' 30113 


Qi = 


9.1186 2.8350" 
2.8350 0.8814. 


,Q 2 = 


' 0.0424 -0.6495" 
-0.6495 9.9576 _ 




Separate opt. 


C X1 = 7.7900, C X1 = 7.7900 

C X2 = 9.9709, C X2 = 9.9686e l0 0067 


Qi = 


7.7900 0" 

,Q 2 = 

oj 


'9.9697 0.03361 
.0.0336 O.OOI2J 


Joint opt. 


C X1 = 9.9905, C X1 = 9.9880 

C X2 = 9.3757, C X2 = 9.3753e- 10 - 4189 


Qi = 


9.9893 " 
0.0013. 


,Q 2 = 


' 8.9702 -1.9067" 
-1.9067 0.4055 _ 





TABLE II: Sum-Rate comparison 





WMMSE 


Separate opt. 


Joint opt. 


Rate-pair 


(2.9220, 1.8401) 


(3.3177,2.0893) 


(3.3307,2.2907) 


Sum-rate 


4.7621 


5.4070 


5.6214 


Improvement 




13.54% 


18.04% 



widely linear precoding, which efficiently maps the proper Gaussian information-bearing signals to the 
improper Gaussian transmitted signals with a given pair of covariance and pseudo-covariance matrices. 
Furthermore, for the two-user SISO-IC, we formulated the optimization problem to characterize the Pareto 
boundary of the achievable rate region via the rate-profile method. Both joint and separate covariance and 
pseudo-covariance optimization algorithms were proposed, both of which outperform the conventional 
proper Gaussian signaling and provide advantages over existing improper Gaussian signaling schemes. 

Appendix A 
Proof of Theorem|4] 

For notational convenience, in this appendix, we use X\ and X2 to represent C Xl and C X2 , respectively. 
First, the following proposition shows that to solve (PI. 7), we may consider exterior solutions only, i.e., 
the solutions at which at least one of the inequality constraints is active. 

Proposition 1. If {X\,X2} is feasible to (PI.7) with \Xi\ < C* and \X2\ < C* 2 , then there exists 
another feasible solution {X[,X 2 } with \X[\ = C* or \X' 2 \ = C* 2 . 

Proof: Let r = min {txtj p^f }• Then r > 1. Define X[ = tX\, X' 2 = tX^. Then the constraints 
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in (S3 1 and (54) are satisfied by X[ and X 2 , i.e., \X[\ < C* and \X' 2 \ < C* 2 . Furthermore, at least one 



of them is satisfied with equality. The constraint in ( |5T| ) is also satisfied since 

ai\h 2 n X[ + h 2 2 X' 2 \ 2 + b x = r^il/ifiXi + h\ 2 X 2 \ 2 + h 



(a) (b) 

< r^axl^Xi + h 2 12 X 2 \ 2 + 61) < r 2 |X 2 | 2 



1 v' 1 2 

l A 2l > 



where (a) comes from r > 1 and 61 > 0, (b) is true since {X\, X 2 } is feasible to (PI. 7). Similarly, 
is also satisfied. Therefore, {X[,X' 2 } is a feasible solution to (P1.7) with at least one of the inequality 
constraints being active. ■ 
Next, we derive Theorem [4] using the Karush-Kuhn-Tucker (KKT) conditions, which are necessary 



optimality conditions for the constrained optimization problem (PI. 7) |52|. For notational convenience, 
denote the inequality constraints ([5T]>— ([54|> by fi < 0, f 2 < 0, h\ < and h 2 < 0, respectively. Denote 
Ai, A2, Hii ^2 as the corresponding dual variables, respectively. Then the Lagrangian function of (P1.7) 
is given by 

L(X 1 ,X 2 ,\ 1 ,\ 2 ,v 1 , f i 2 ) = \ 1 {a 1 (\h 11 \ 4 X*X 1 + \h 12 \ 4 X*X 2 
+ 2^{h 2 \h 2 2 XlX 2 }) + 61 - X* 2 X 2 } + \ 2 {a 2 {\h 2 i\ A XlX l 
+ \h 22 \ 4 X* 2 X 2 + 2^{hf l h 2 22 XlX 2 }) + b 2 - XlXi} 



+ nx(xzx x - c: 2 ) + ii 2 {x* 2 x 2 - a 



*2\ 
x 2 ) 



(57) 



For {X\,X 2 } to be a solution to (PI. 7), the following KKT conditions must be satisfied: 
1) Dual feasibility: Ai > 0, A 2 > 0, m > 0, fx 2 > 0. 



2) Zero derivative: The derivatives of the Lagrangian function ( |57j ) with respect to the primal variables 

=> - X 2 [(ai\h 12 \ 4 - l)Ai + A 2 a 2 |/i 22 | 4 + fJa] 



are zero [55 1 



dL 



dX* 



dL 

dX* 



= Xi(Ai a 1 h 1 * 2 h 11 +A 2 a 2 h 22 h 21 ) 
Vi v 2 

=^Xi [(a 2 |/i 2 i| 4 - 1)A 2 + Aiai|/in| 4 + m] 



Cl 



-X 2 (Xi aih^h l2 +A 2 a 2 h 2 \h 22 ) 
V{ ~V> 



-X 2 c 2 =X 1 (X 1 V 1 + X 2 V 2 ), 
X x c x = -X 2 (\iVi + X 2 V 2 *) 



(58) 
(59) 
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3) Complementary slackness: \\f\ = 0,A2/2 = 0,//i/ti = 0,^2^2 = 0- 
As discussed previously, without loss of generality, X\ can be assumed to be a nonnegative real number. 
The following cases are then considered to derive the possible phases of X2'. 

• Case I: /1 = and /2 / 0. Then from the complementary slackness condition, we have Ai > and 



A2 = 0. Substituting them into (59 1, we have 



Xi(Aiai|/in| 4 + yLi) 

X2 = w — Vl - 

Since Ai > 0, > 0, a\ > and X\ > 0, the phase 9 of X2 equals to that of V\ rotated by it, 
which is 7T + 2(0n — 012) since V\ = aihf^hfi- 



Case II: /1 / and f 2 = 0. Then Ai = and A2 > 0. Similarly, by using ( [58] ), we have 

9 = 7T + 2(021 - 02 2 )- 



Case III: fi^O and f 2 / 0, then Ai = and A 2 = 0. By substituting them into ([58]) and ([59]), we 
have [12X2 = and [i\X\ = 0. This means either X2 = 0, X\ = 0, or [i\ = 0, fi 2 = 0. Furthermore, 
Proposition [T] suggests that we may consider the exterior solutions only, i.e., either hi = or /12 = 
is satisfied. Thus, /xi > or fi 2 > can be assumed. Therefore, case III can be ignored without 
loss of optimality. 

Case IV: f\ = and f 2 = 0, then Ai > and A2 > 0. In this case, 9 belongs to the solution set for 
the equations given in ( [55] ) and ([56]), which are obtained by satisfying the constraints in ( [51] ) and 



( [52] ) with equality. p5\ and ( [56] ) correspond to \X±\ = C* and \X2\ = C* 2 , respectively, which can 
be assumed without loss of generality due to Proposition [T] 
This completes the proof of Theorem |4] 

Appendix B 
Solving O.4 and © b in Theorem[4] 

In this appendix, we show the steps to solve 0.4. ©0 can be obtained similarly. The unknown variables 
in ( |55] l are 9 and t. After some manipulations, ([55]) can be written as 



tcosrj + dit 2 + d 2 = (60) 
t cos(rj + u) + d 3 t 2 + d± = (61) 



where 

A J. J. J. J. \ A 



UJ = 2(022 + 011 - 012 - 021 ), rj = 9 + 2(012 ~ 01l) (62) 



dl — ~ 1 r ' i'i ) ^2 



ai|/ii 2 | 4 -l A ai|/in| 4 C* 2 + &i 



2ai|/in| 2 |/ii 2 | 2 Q i ' 2ai|/in| 2 |/ii 2 | 2 a 
12 f„„i/,„,|4_ nr-*2 



4 |fr 22 | 2 (a 2 |/i 2 i| 4 -l)Qf + 6 2 

3 2]fc 21 ] 2 C*' 4 2a2|/ l2 i| 2 |/i22| 2 Ci i ' 
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From doTh, we have 



t sin 77 sin u; = t cos 77 cos u; + d%t + (I4 ^ 

t 2 (l — cos 2 77) sin 2 u) = (t cos 77 cos a; + d^t 2 + d^) 2 



(63) 



Solving cos 77 from ( |60| ), we have 

cos 77 = - (c?it 2 + d 2 ) /t. (64) 
Substituting (64 ) into ([63]> gives the following fourth order polynomial equation with respect to t: 



[t 2 — (d\t 2 + d2) 2 ] sin 2 oj = [(d% — d\ cosui)t 2 + d& — d2 cos uj} 2 . 

Since the above equation only has t 2 terms, it can be transformed to the following quadratic equation by 
setting z = t 2 , 

eiz 2 + e 2 z + e 3 = 0, (65) 

where e\ =d\ + d\ — 2d\d% cos uj, 

e 2 =2{d\d2 + d^di) — 2(d\d4 + d 2 d^) cos 00 — sin 2 uj, 

£3 + d 2 — 2d 2 d^ cos uj. 

Then z can be easily solved. Since z = t 2 and t is the magnitude of C X2 , only the solutions of z that are 
real and satisfy < z < C* 2 need to be kept, whereby the values for t are obtained. For those values of t 
satisfying | (dit 2 + d2)/t\ < 1, we can get the value for 77 based on (64l, i.e., 77 = arccos[— (d\t 2 + d2)/t\ 
or 77 = 2tt — arccos[— (d\t 2 + d 2 )/t}. Then 9 can be obtained from (62 1. If no such solutions exist, then 
0,4 is set to empty. 

Appendix C 
An Upper Bound R v for (PI. 6) 

In this section, we derive an upper bound Rjj for the optimal solution R* to the pseudo-covariance 
optimization problem (PI. 6). Based on our discussions in Section |Vj the target rate R should be feasible 



to (PI. 7). From the constraints in (51 1 and (54i, it can be seen that R should satisfy: 



61 < a, r < ct 



±2 



(66) 



where b 



(l-l/flQC? 

I^12| 4 



and fii = < 
1 



2on(R r*)_ gy substituting b\ into ( [66] ), we have 

h^C* 2 a A + 2a 2 \h l2 \ 2 Ct 



h ~ 1 ci 



x 2 



where for the last equality, we have used the identity C Sl 

>12 



\h\2\ 2 C* +a 2 . Then we have 



&<1 + 



^12| 4 C* 2 



a 2 (^ + 2\h 12 \ 2 C* 2 ) 



(67) 
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By substituting f3\ into (67i, we obtain 



1 



|/ll2l 4 C* 2 



x 2 



2«i b V o- 2 (o- 2 + 2\h 12 \ 2 C* 2/ 



(68) 



A similar upper bound can be obtained based on the constraints in p2\ and ( |53j ). Thus the upper bound 
for R* in (PI. 6) is given by 



Ru = r* + min 



U. 

\2ai 
1 



log 1 + 



|/il2| 4 C^ 



a 2 (a 2 + 2|/ ll2 |2Q, 



, A l fe 2l| 4 C* 2 

2a 2 g V i + o- 2 ^ 2 + 2|/ l21 |2Q i ; 



(69) 



References 

[1] C. E. Shannon, "Two-way communication channels," in Proc. 4th Berkeley Symp. Math. Stat, and Prob., Berkeley, CA, 
1961, pp. 611-644. 

[2] A. B. Carleial, "A case where interference does not reduce capacity," IEEE Trans. Inf. Theory, vol. 21, no. 5, pp. 569-570, 
Sep. 1975. 

[3] T. S. Han and K. Kobayashi, "A new achievable rate region for the interference channel," IEEE Trans. Inf. Theory, vol. 27, 
no. 1, pp. 49-60, Jan. 1981. 

[4] R. Etkin, D. Tse, and H. Wang, "Gaussian interference channel capacity to within one bit," IEEE Trans. Inf. Theory, 

vol. 54, no. 1, pp. 5534-5562, Dec. 2008. 
[5] X. Shang, G. Kramer, and B. Chen, "A new outer bound and noisy-interference sum-rate capacity for the Gaussian 

interference channels," IEEE Trans. Inf. Theory, vol. 55, no. 2, pp. 689-699, Feb. 2009. 
[6] A. S. Motahari and A. K. Khandani, "Capacity bounds for the Gaussian interference channel," IEEE Trans. Inf. Theory, 

vol. 55, no. 2, pp. 620-643, Feb. 2009. 
[7] M. Chiang, R Hande, T. Lan, and C. W. Tan, Power Control in Wireless Cellular Networks. Foundations and Trends in 

Networking, 2008. 

[8] A. Gjendemsj0, D. Gesbert, G. E. 0ien, and S. G. Kiani, "Binary power control for sum rate maximization over multiple 

interfering links," IEEE Trans. Wireless Commun., vol. 7, no. 8, pp. 3164-3173, Aug. 2008. 
[9] Z. Q. Luo and S. Zhang, "Dynamic spectrum management: complexity and duality," IEEE J. Set Topics Signal Process., 

vol. 2, no. 1, pp. 57-73, Feb. 2008. 
[10] F. R. Farrokhi, K. J. R. Liu, and L. Tassiulas, "Transmit beamforming and power control for cellular wireless systems," 

IEEE J. Sel. Areas Commun., vol. 16, no. 8, pp. 1437-1450, Oct. 1998. 
[11] E. Visotsky and U. Madhow, "Optimal beamforming using transmit antenna arrays," in Proc. IEEE Veh. Technol. Conf., 

vol. 1, Houston, Texas, May 1999, pp. 851-856. 
[12] M. Schubert and H. Boche, "Solution of the multiuser downlink beam-forming problem with individual SINR constraints," 

IEEE Trans. Veh. Technol, vol. 53, no. 1, pp. 18-28, Jan. 2004. 
[13] H. Dahrouj and W. Yu, "Coordinated beamforming for the multicell multi-antenna wireless system," IEEE Trans. Wireless 

Commun., vol. 9, no. 5, pp. 1748-1759, May 2010. 
[14] B. Song, R. Cruz, and B. Rao, "Network duality for multiuser MIMO beamforming networks and applications," IEEE 

Trans. Commun., vol. 55, no. 3, pp. 618-630, Mar. 2007. 
[15] A. Wiesel, Y. C. Eldar, and S. Shamai (Shitz), "Linear precoding via conic optimization for fixed MIMO receivers," IEEE 

Trans. Signal Process., vol. 54, no. 1, pp. 161-176, Jan. 2006. 



July 24, 2012 



DRAFT 



29 



[16] M. Bengtsson and B. Ottersten, "Optimal downlink beamforming using semidefinite optimization," in Proc. 37th Allerton 

Conf. on Commun., Control, and Computing, Mar. 1999, pp. 987-996. 
[17] Y. F. Liu, Y. H. Dai, and Z. Q. Luo, "Coordinated beamforming for MISO interference channel: complexity analysis and 

efficient algorithms," IEEE Trans. Signal Process., vol. 59, no. 3, pp. 1142-1156, Mar. 2011. 
[18] J. Huang, R. A. Berry, and M. L. Honig, "Distributed interference compensation for wireless networks," IEEE J. Sel. Areas 

Commun., vol. 24, no. 5, pp. 1074-1084, May 2006. 
[19] R. Zakhour and D. Gesbert, "Coordination on the MISO interference channel using the virtual SINR framework," in Proc. 

1TG/IEEE Work-shop Smart Antennas, Feb. 2009. 
[20] S. Ye and R. S. Blum, "Optimized signaling for MIMO interference systems with feedback," IEEE Trans. Signal Process., 

vol. 51, no. 11, pp. 2839-2848, Nov. 2003. 
[21] H. Sung, K. J. Lee, S. H. Park, and I. Lee, "Linear precoder designs for Jf-user interference channels," IEEE Trans. 

Wireless Commun., pp. 291-301, Jan. 2010. 
[22] L. R Qian, Y. J. Zhang, and J. Huang, "MAPEL: achieving global optimality for a non-convex wireless power control 

problem," IEEE Trans. Wireless Commun., vol. 8, no. 3, pp. 1553-1563, Mar. 2009. 
[23] E. A. Jorswieck and E. G. Larsson, "Monotonic optimization framework for the two-user miso interference channel," IEEE 

Trans. Commun., vol. 58, no. 7, pp. 2159-2168, Jul. 2010. 
[24] L. Liu, R. Zhang, and K. C. Chua, "Achieving global optimality for weighted sum-rate maximization in the A"-user 

Gaussian interference channel with multiple antennas," IEEE Trans. Wireless Commun., vol. 11, no. 5, pp. 1933-1945, 

May 2012. 

[25] W. Utschick and J. Brehmer, "Monotonic optimization framework for coordinated beamforming in multicell networks," 

IEEE Trans. Signal Process., vol. 60, no. 4, pp. 1899-1909, Apr. 2012. 
[26] E. Bjornson, G. Zheng, M. Bengtsson, and B. Ottersten, "Robust monotonic optimization framework for multicell MISO 

systems," IEEE Trans. Signal Process., vol. 60, no. 5, pp. 2508-2523, May 2012. 
[27] S. S. Christensen, R. Agarwal, E. D. Carvalho, and J. M. Cioffi, "Weighted sum-rate maximization using weighted MMSE 

for MIMO-BC beamforming design," IEEE Trans. Wireless Commun., vol. 7, no. 12, pp. 4792-4799, Dec. 2008. 
[28] M. Razaviyayn, M. Sanjabi, and Z. Q. Luo, "Linear transceiver design for interference alignment: complexity and 

computation," IEEE Trans. Inf. Theory, vol. 58, no. 5, pp. 2896-2910, May 2012. 
[29] G. Scutari, P. Palomar, and S. Barbarossa, "Competitive design of multiuser MIMO systems based on game theory: a 

unified view," IEEE J. Sel. Areas Commun., vol. 26, no. 7, pp. 1089-1103, Sep. 2008. 
[30] A. Host-Madsen and A. Nosratinia, "The multiplexing gain of wireless networks," Int. Symp. on Inf. Theory and Its 

Applications, pp. 2065-2069, 4-9 Sept. 2005. 
[31] V. R. Cadambe and S. A. Jafar, "Interference alignment and degrees of freedom of the if -user interference channel," IEEE 

Trans. Inf. Theory, vol. 54, no. 8, pp. 3425-3441, Aug. 2008. 
[32] C. M. Yetis, T. Gou, S. A. Jafar, and A. H. Kayran, "On feasibility of interference alignment in MIMO interference 

networks," IEEE Trans. Signal Process., pp. 4771 - 4782, Sep. 2010. 
[33] J. Park, Y. Sung, and H. V. Poor, "On beamformer design for multiuser MIMO interference channels," 2010, [Online]. 

Available: http://arxiv.org/abs/101 1.6121. 
[34] S. W. Peters and R. W. Heath Jr., "Cooperative algorithms for MIMO interference channels," IEEE Trans. Veh. Technol, 

vol. 60, no. 1, pp. 206 - 218, Jan. 2011. 
[35] D. A. Schmidt, C. Shi, R. A. Berry, M. L. Honig, and W. Utschick, "Minimum mean squared error interference alignment," 

IEEE Asilomar Conference on Signals, Systems and Computers (ACSSC), pp. 1106 - 1110, Nov. 2009. 
[36] K. Gomadam, V. R. Cadambe, and S. A. Jafar, "A distributed numerical approach to interference alignment and applications 

to wireless interference networks," IEEE Trans. Inf. Theory, vol. 57, no. 6, pp. 3309-3322, May 2011. 



July 24, 2012 



DRAFT 



30 



[37] M. Charafeddine, A. Sezgin, and A. Paulraj, "Rates region frontiers for n-user interference channel with interference as 

noise," in Proc. Allerton Conference, Sep. 2007. 
[38] E. Jorswieck, E. Larsson, and D. Danev, "Complete characterization of the Pareto boundary for the MISO interference 

channel," IEEE Trans. Signal Process., vol. 56, no. 10, pp. 5292-5296, Oct. 2008. 
[39] R. Zhang and S. Cui, "Cooperative interference management with MISO beamforming," IEEE Trans. Signal Process., 

vol. 58, no. 10, pp. 5450-5458, Oct. 2010. 
[40] X. Shang, B. Chen, and H. V. Poor, "Multiuser MISO interference channles with single-user detection: optimality of 

beamforming and the achievable rate region," IEEE Trans. Inf. Theory, vol. 57, no. 7, pp. 4255 - 4273, Jul. 2011. 
[41] R. Mochaourab and E. A. Jorswieck, "Optimal beamforming in interference networks with perfect local channel 

information," IEEE Trans. Signal Process., vol. 59, no. 3, pp. 1128-1141, Mar. 2011. 
[42] F. D. Neeser and J. L. Massey, "Proper complex random processes with applications to information theory," IEEE Trans. 

Inf. Theory, vol. 39, no. 4, pp. 1293 - 1302, Jul. 1993. 
[43] P. J. Schreier and L. L. Scharf, Statistical Signal Processing of Complex-Valued Data: The Theory of Improper and 

Noncircular Signals. Cambridge (UK): Cambridge Univ. Press, 2010. 
[44] G. Taubock, "Complex-valued random vectors and channels: entropy, divergence, and capacity," IEEE Trans. Inf. Theory, 

vol. 58, no. 5, pp. 2729-2744, May 2012. 
[45] S. A. Jafar, Interference Alignment: A New Look at Signal Dimensions in a Communication Network. Foundations and 

Trends in Communications and Information Theory, 2010, vol. 7, no. 1. 
[46] V. R. Cadambe, S. A. Jafar, and C. Wang, "Interference alignment with asymmetric complex signaling - settling the 

Host-Madsen-Nosratinia conjecture," IEEE Trans. Inf. Theory, pp. 4552 - 4565, Sep. 2010. 
[47] Z. K. M. Ho and E. Jorswieck, "Improper Gaussian signaling on the two-user SISO interference channel," in Proc. IEEE 

ISWCS, 2011. 

[48] S. H. Park, H. Park, and I. Lee, "Coordinated SINR balancing techniques for multi-cell downlink transmission," in Proc. 
VTC 2010-fall, 2010. 

[49] P. J. Schreier and L. L. Scharf, "Second-order analysis of improper complex random vectors and processes," IEEE Trans. 

Signal Process., vol. 51, no. 3, pp. 714-725, Mar. 2003. 
[50] B. Picinbono and P. Chevalier, "Widely linear estimation with complex data," IEEE Trans. Signal Process., vol. 43, no. 8, 

pp. 2030-2033, Aug. 1995. 

[51] G. Gelli, L. Paura, and A. R. P. Ragozini, "Blind widely linear multiuser detection," IEEE Commun. Lett., vol. 4, no. 6, 
pp. 187-189, Jun. 2000. 

[52] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge, U.K.: Cambridge Univ. Press, 2004. 

[53] R. L. Streit and A. H. Nuttall, "A note on the semi-infinite programming approach to complex approximation," Math. 

Computat., vol. 40, no. 162, pp. 599-605, Apr. 1983. 
[54] M. Grant and S. Boyd, CVX: Matlab software for disciplined convex programming, version 1.21, http://cvxr.com/cvx. 
[55] A. Hj0rungnes, Complex-Valued Matrix Derivatives With Applications in Signal Processing and Communications. 

Cambridge (UK): Cambridge Univ. Press, 2011. 



July 24, 2012 



DRAFT 



